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■ Abstract 

o: 

We investigate the resonant leptogenesis scenario in the minimally U{1)b~l extended 
standard model with minimal flavor violation. In our model, the U{1)b-l gauge symmetry is 
broken at the TeV scale and standard model singlet neutrinos gain Majorana masses of order 
TeV. In addition, we introduce a flavor symmetry on the singlet neutrinos at a scale higher 
5^ ! than TeV. The flavor symmetry is explicitly broken by the neutrino Dirac Yukawa coupling, 
which induces splittings in the singlet neutrino Majorana masses at lower scales through renor- 
malization group evolutions. We call this setup "minimal flavor violation" . The mass-splittings 
are proportional to the tiny Dirac Yukawa coupling, and hence they automatically enhance the 
CP asymmetry parameter necessary for the resonant leptogenesis mechanism. In this paper, 
we calculate the baryon number yield by solving the Boltzmann equations, including the effects 
of U{1)b-l gauge boson that also has TeV scale mass and causes washing-out of the singlet 
neutrinos in the course of thermal leptogenesis. The Dirac Yukawa coupling for neutrinos is 
fixed in terms of neutrino oscillation data and an arbitrary 3x3 complex-valued orthogonal 
matrix. We show that the right amount of baryon number asymmetry can be achieved through 
thermal leptogenesis in the context of the minimal flavor violation with singlet neutrinos and 
U{1)b-l gauge boson at the TeV scale. These particles can be discovered at the LHC in the 
near future. 



1 Introduction 



The baryon number asymmetry of the Universe has been measured to be [T] 

— = (0.88 ±0.02) X 10^1° . (1) 
s 

Since the electroweak baryogenesis scenario based on the standard model (SM) is incapable of 
explaining this amount, we need an extension of SM. Leptogenesis [2] is an appealing scenario, 
as it explains the tiny neutrino masses and the baryon number asymmetry simultaneously 
with heavy singlet neutrinos [3]. What is problematic in the leptogenesis scenario is that it 
usually cannot be tested experimentally; the efficiency of leptogenesis is proportional to CP 
violation which appears in interference terms between tree and one-loop diagrams involving 
singlet neutrinos. Therefore sufficiently large neutrino Yukawa coupling is required in order 
that the efficiency be high enough to produce the observed amount of baryon asymmetry. This 
in turn implies that the singlet neutrinos must be heavy enough to explain the tiny neutrino 
masses, and hence they are beyond the energy scale of collider experiments. 

However, if singlet neutrinos have nearly degenerate Majorana masses, resonance in loop di- 
agrams enhances the efficiency of leptogenesis, making it possible to produce the right amount 
of baryon asymmetry even with small neutrino Yukawa couplings, which is called "resonant 
leptogenesis" With the resonant leptogenesis mechanism, models with TeV scale singlet 
neutrinos can explain both the baryon number asymmetry and the light neutrino masses. Still, 
collider experiments cannot test such models because the singlet neutrinos couple to SM parti- 
cles only through the small Yukawa coupling, and hence they cannot be produced at colliders. 
If there is a new gauge interaction that is broken at TeV scale and the singlet neutrinos are 
charged under it, we may be able to produce the TeV scale singlet neutrinos at colliders and 
even measure the amount of CP violation through their decays j5]. Another problem of mod- 
els with TeV scale singlet neutrinos is its naturalness. Since successful resonant leptogenesis 
requires highly degenerate Majorana masses, we have to explain the origin of the degeneracy 
from some new symmetry or mechanism. 

In this paper, we propose a model that solves these problems and show that it can ex- 
plain both the light neutrino masses and the baryon number asymmetry without losing its 
experimental accessibility and without requiring unnatural fine-tuning. In addition to the SM 
content, our model has U{1)b~l gauge symmetry and contains three singlet neutrinos which 
are necessary for its anomaly cancellation. This gauge symmetry is broken at the TeV scale so 
that U{1)b~l gauge boson and the singlet neutrinos obtain masses of order TeV [6]. We further 
introduce a flavor symmetry on the singlet neutrinos at an energy scale above TeV, which is 
broken explicitly by the neutrino Dirac Yukawa coupling. 
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Because of the flavor symmetry of the singlet neutrinos, the differences in their Majorana 
masses arise solely from the neutrino Dirac Yukawa coupling through renormalization group 
(RG) evolutions, which we call "minimal flavor violation" for singlet neutrinos. The ratios of 
the mass differences over the mass scale are therefore of the same order as the neutrino Dirac 
Yukawa coupling. It turns out that such mass difference is of an appropriate order to produce 
the right amount of the baryon number asymmetry with resonant enhancement, and no fine- 
tuning is required in our setup. Moreover, since the differences in the Majorana masses are 
related to the neutrino Dirac Yukawa coupling, our model is more restrictive than the original 
resonant leptogenesis scenario, where the mass difference is a free parameter. 

Two of the authors have discussed in ref. [7] resonant leptogenesis in the presence of U{1)b-l 
gauge symmetry that is broken at TeV scale. In that paper, the tiny difference in the Majorana 
masses of two singlet neutrinos is put in by hand, whereas in our paper it arises from RG 
evolutions involving the neutrino Dirac Yukawa coupling. 

Resonant leptogenesis with a flavor symmetry of singlet neutrinos at a high scale has been 
first proposed in ref. [5]. However the models in that paper does not contain U{1)b~l gauge 
symmetry that is broken at TeV scale because its motivation lies in lowering the reheating 
temprature of the Universe rather than observing singlet neutrinos by collider experiments. We 
expect that the washing-out of singlet neutrinos via U{1)b~l gauge interaction in the course 
of thermal leptogenesis reduces the resultant baryon number asymmetry, making it difficult 
to produce the right amount of baryon asymmetry. This is particularly so when the mass of 
U{1)b~l gauge boson is larger than twice the Majorana masses, because in this case on-shell 
U{1)b-l gauge boson can mediate the interactions between singlet neutrinos and SM particles. 
From the point of view of collider experiments, however, this case is favorable because the 
singlet neutrinos can be produced by the decay of U{1)b^l gauge boson. We will see that even 
with the U{1)b-l gauge interaction at the TeV scale, we can obtain enough baryon asymmetry 
through resonant leptogenesis with a flavor symmetry of singlet neutrinos at a high scale. 

Another feature of our study is that we adopt Casas-Ibarra parametrization [9] for the 
neutrino Dirac Yukawa coupling. When the charged lepton Yukawa coupling is neglected, the 
processes of leptogenesis depend only on (Yj)ai(Y^)ib, where Y^, denotes the neutrino Dirac 
Yukawa coupling, i the flavor index of lepton doublets and a, b that of singlet neutrinos. In 
Casas-Ibarra parametrization, this quantity is parametrized in terms of the light neutrino mass 
eigenvalues and an arbitrary complex-valued orthogonal matrix, independently of the neutrino 
flavor-mixing matrix. This parametrization thus makes it clear whether the right amount of 
baryon number asymmetry can be produced in a given setup. 

In Section 2, we describe the model which contains U{1)b-l gauge symmetry and a flavor 
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Field 


SUi2)L 


f/(l)y 


U{1)b--l 


U 


2 


-1/2 


-1 


Na 


1 





-1 


H 


2 


-1/2 





$ 


1 





+2 



Table 1: Field content in the minimally U{1)b-l extended standard model. 

symmetry of singlet neutrinos. In Section 3, we discuss its RG evolutions, which are essential 
in our model as they give rise to the small mass-splittings of Majorana masses necessary for 
resonant leptogenesis. In Section 4, we review on the Boltzmann equations that calculate the 
baryon number yield from leptogenesis. In Section 5, we calculate the baryon number yield for 
various values of singlet neutrino masses and U{1)b~l gauge boson mass, and study whether 
the observed amount of the baryon asymmetry can be produced. Section 6 is devoted to the 
conclusions and discussions. 



2 Setup 

We consider the minimally U{1)b~l extended standard model, where the gauge symmetry is 
U{1)y X U{1)b~l X SU{2)l x SU{?>)c and right-handed singlet neutrinos, Na {a = 1,2,3), and 
a scalar field, $, are added to the SM matter content. $ obtains VEV and breaks the U{1)b-l 
gauge symmetry through the Higgs mechanism, giving rise to the U{1)b~l gauge boson mass 
and the Majorana mass term for the singlet neutrinos. We denote this VEV by vb-l- We 
are especially interested in the case where vb-l is at the TeV scale, and hence the U{1)b-l 
gauge boson and the singlet neutrinos are observable at collider experiments. In this case, the 
components of Yi, are smaller than about 10~^. The field content is summarized in Table 1. We 
assume that the physical scalar particle that originates from $ after the U{1)b-l symmetry 
breaking is heavier than the singlet neutrinos, and we may neglect it in the discussion of thermal 
leptogenesis. 

We take the combination of the f/(l)y and U{1)b~l gauge bosons such that the mixing of 
their kinetic terms are zero throughout the renormalization group (RG) evolutions. Then the 
convariant derivative, Z)^, can be written as 

D, = d, ~ ^^9lYB, - i{9Y + gB-LYB-L)B'^ - 192^^ " ^^?3T°G^ , (2) 
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where the couphng gY reflects the mixing of the U{1)y and U{1)b~l gauge bosons. Q^Qb-l 
are defined so that we have ^ = at the U{1)b-l breaking scale. Since $ has B — L charge 2, 
its VEV induces the mass of the U{1)b~l gauge boson as 

Mz> = 2V2gB^LVB~L . (3) 
The Dirac and Majorana Yukawa couplings for the neutrinos are given by 

C 3 ~{Yu)^a HL\Na - h.C. - ]^{YM)ab ^N,N, - h.C. (4) 

Note that Ym is a symmetric matrix with complex values. The term \{YM)ah ^NaNf, becomes 
the Majorana mass term, 

^VB-L{YM)ab NaNb , (5) 

when $ obtains VEV. 

The light neutrino mass matrix takes the form: 



[m 



UjlJ 



iY,)^a (yj),, (6) 



vb-l 



We move to the flavor basis where the charged leptons and the singlet neutrinos are in their 
mass eigenstates. On this basis, we have 

/ 

mlm^ = Umns | Ul^^g , (7) 

\ mj 

2/1 
Ym = I 2/2 
2/3 

where (mj)^'s denote the eigenvalues of the light neutrino mass squared, Umns denotes the 
neutrino flavor-mixing matrix and ?/j's denote the eigenvalues of Ym- We see from (6), (7) and 
(8) that Yi, on this basis can generally be written as 

\ /v^r \ 

iY,)^a = Ul,^s (m2)V4 Q O , (9) 

V {mlf/' I V ^ 

where O denotes an arbitrary orthogonal matrix with complex values, which appears in Casas- 
Ibarra parametrization [9] . We hereafter call O "Yukawa kernel matrix" . 
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Let us introduce a symmetry of two generations of singlet neutrinos above an energy scale 
A, which is assumed higher than the TeV scale. Then we may write 



in the basis where Ym is diagonal. This is the key assumption that realizes the resonant 
leptogenesis without fine-tuning of the parameters. Below the scale A, Yi, explicitly breaks the 
flavor symmetry and comes into effect through RG evolutions. The difference between the first 
two components of Ym at lower scales is proportional to YjYi^, and hence the singlet neutrinos 
obtain Majorana masses whose difference between the first two components is proportional to 
YjYy when $ gains VEV. This mass difference is essential for successful resonant leptogenesis. 

Our flavor symmetry is inspired by a class of A4 symmetric models where the three flavors 
of the singlet neutrinos transform as 3 representation of A4 group (see e.g. Table 1 of ref. [10] 
for a list of such models), although we do not specify the details of the flavor symmetry. If the 
singlet neutrinos are in 3 representation and the B — L breaking field $ is in 1 representation 
of ^4 group, the Majorana Yukawa coupling will be common for the three flavors. In realistic 
/I4 models, especially when one takes into account the observed value of the neutrino mixing 
angle 613 [H], A4 symmetry should be spontaneously broken below some scale. In our model, 
the spontaneous breaking of a flavor symmetry is effectively described by the Dirac Yukawa 
coupling that explicitly violates the flavor symmetry and that comes into play below the scale 
A. 



In the following discussion, we set the U{1)b-l breaking scale, Mq, at \y\vB~L, where y is 
of the order of the eigenvalues of Ym- The results do not change with 0(1) variation of the 
U{1)b-l breaking scale. We further make an approximation that the U{1)b-l gauge boson 



The RG equations for the U(l) gauge couplings, gi,gB-L,g, are given as follows, where Mq 




(10) 



3 RG Evolutions 



mass, Mz' = 2\/2gB-LVB-L, is of the same order as Mq. 
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is the U{1)b^l breaking scale: 
2 d 53 o 

2 d 41 o „ , , 

IGtt /i-j-fi-i = —5-1 lor /i > mf , (11) 

dyU 1(J 

16vrV^^7B-L = 12f7|_,. + 2 ■ y^?!.,.^ + ^Sb-l^ for > Mq , (12) 

9 d _ 41„,„o 6 2n „ 16 ,^2 3 2x ,^ 9 _ ^ 

167rV:^^7 = ir^?(r + 7^?i) + ^ ■ ^93-^9^ + -gl) + 1^9b-l9 for /i > Mq. 
d/x D 5 3 5 

(13) 

The RG equations for g2igz are the same as in the standard model. 

The RG equations for the Yukawa couplings Ym, Y^Y^, and the top Yukawa coupling, ?/f , are 
given as : 

167rV:^(FM)a6 = {YX)ac{YM)cb + {Ym) aciY^u) cb 



d/i 



+ {YMY;,YM)ab + ]^iYM)abtT[YMYl] - Q{g'^fiYM)ab for /i > Mq 



(14) 



lQ'K^^^^^{YX)al> = {YX)ab ( - - ^gj - Ugl^^ - - 6gB-L~g 



+ ^{YMY*jYX)ah + ^{YXYMY*,i)ab for /i > Mq , (15) 



d 9 . 17 9 

16vrV^yt = 2^** ~ ^* ( + + 8^3 ) for mj < ^ < Mq 

d 9 . 17 9 17 2 5_ 

167rV^I/t = ^Vt - ( ^gf + |£/2 + %l + -j^f + 3^Ll + ^ggs^L ) for /i > Mq. 

(16) 

In the RG equations above, we neglect the Yukawa couplings except the top Yukawa coupling. 
Note that, when the charged lepton Yukawa couplings are neglected, the RG evolution for Yi, 
can be expressed in terms of that for YjY^. 
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The boundary conditions for the RG equations are 
/ y \ 

Ym[A] = { y \ , (17) 



vb~l 

y2 




'ml 

y}y,[Mq] ~ I I oM o I o 



^Jml 




(18) 

~g[Mo] = 0. (19) 

Note that Y^Y^ is independent of the neutrino flavor-mixing matrix, Umns- In the first and 
second equations, y,yz and y,yz are related through RG evolutions from the scale A to Mq. In 
order to obtain their relation, we solve eq. (14) neglecting the term Y^Y^. 

We numerically solve the RG equations using the following best fit values (with 3cr ranges) 
from the Particle Data Group [12], which were derived assuming non-zero 6'i3: 

ml -ml = 7.58 tofg x 10"^ eV^ , (20) 
\ml-ml\ = 2.35 +°;| x lO'^ eV^ , (21) 

Our analysis is done for both the normal hierarchy (m| > ml) and the inverted hierarchy 
(m| < ml) of the light neutrino masses. We fix the mass of the lightest neutrino at zero, i.e., 
we set ml = for the normal hierarchy case and mg = for the inverted hierarchy case. The 
Yukawa kernel matrix O can be any 3x3 orthogonal matrix. However, since the purpose of this 
paper is to show that successful leptogenesis can be realized in our framework, we here restrict 
our analysis to a specific choice of the matrix O that can be expressed as 



O = I cos a sin a for the normal hierarchy case. 




O = —sin a cos a for the inverted hierarchy case, (22) 



where a takes a complex value. With these assumptions, the upper- left 2x2 part of yJyi/[Mo] 
takes the following simple form for the normal hierarchy case: 

I . / cos a* -sin a* \ / ^/m^ \ f cos a sina \ \y\vB-L 



sma" cos a" I \ (J Wmi J \ — sma cos a / v'^ 
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For the inverted hierarchy case, y^m| and \fn\^ in the above formula are replaced with y^mf 
and \frn^. 

We calculate the mass spectrum for various values of a and substitute it into the Boltzmann 
equations that calculate the baryon number yield. Since we are assuming a symmetry on only 
two generations of singlet neutrinos, it is natural that the other singlet neutrino has distinc- 
tively different Majorana mass. This singlet neutrino then does not contribute to resonant 
leptogenesis. We here neglect this singlet neutrino and make a two-flavor approximation, i.e., 
we only consider the two generations of singlet neutrinos, on solving the Boltzmann equations. 



4 Baryon Number Asymmetry from Leptogenesis 

At high temperatures, the singlet neutrinos as well as the SM particles are in thermal 
equilibrium. (Remember that the IJ{\)b~l gauge boson mediates interactions between the 
singlet neutrinos and other particles.) As the Universe cools down and the temperature becomes 
below the order of Mq, the singlet neutrinos decouple from the thermal bath and decay into 
the Higgs doublet and the SM leptons. Because of the CP phases of Y^, the decays of the 
singlet neutrinos exhibit a CP asymmetry, which gives rise to non-zero lepton number density, 
or equivalently non-zero B — L charge density. The lepton number density is converted into 
the baryon number asymmetry through the sphaleron process [151 IIS] when the temperature 
of the Universe is higher than the decoupling temprature of the sphaleron process T^ph ~ 150 
GeV [n]. 

To estimate the resultant baryon number asymmetry, we solve the Boltzmann equations 
governing the number densities of the particles concerned [13]. We normalize the temperature, 
T, of the radiation-dominated Universe as 

Z = yVB-L / T , 

where yvs-i approximately equals to the mass of the two singlet neutrinos that are relevant to 
leptogenesis. We respectively denote the yields, n/s {n: number density of a particle species, 
s: entropy density), of the singlet neutrinos A''^ (a = 1, 2) and the B — L charge by Fat^, Yb-l- 
The CP asymmetry parameter for the decay of singlet neutrino Na is defined as 

T.,[nNa^L,W)-nN,^L]H)] 

E,[r(A^a ^ L,H^) + r(A^„ ^ l]h)] ' 
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H 



H 



H 



Figure 1: Decay of a singlet neutrino at tree and one- loop levels. 



Non-zero ea comes from the interference of the diagrams in Figure 1. It is generally given by 



E 



Mb Mb 



s, 



(25) 



bb 



where Mq, Fq denote the Majorana mass and the decay width of the singlet neutrino Na and 
we have F^ = Ma{Y^Yy)aa/^T^ ■ H reflects the interference of the left and the middle diagrams 
in Figure 1 and that of the left and the right diagrams. They are given by 



Sb 



'Ml 



1 + 



Ml 
Ml 
-M^ 



log 1 + 



Ml 
Ml 



(26) 
(27) 



{Ml - MIY + MlTl ■ 

In our model, the mass difference of the relevant two singlet neutrinos, M| — , comes from 
the RG evolutions from the scale A down to Mq. Therefore it is estimated as 



Ml - Ml oc YX {yvB^L) 



Since we are assuming i/vb^l ~TeV, it follows from eq. (9) that the components of Y^, is at least 
as small as (Am^^^/TeV)^/'^ ~ 10~^. Then 5*1, 5*2 are enhanced by the factor of (l/YjY^) ~ 10*, 
whereas Vi, Vi remain 0(1). In eq. (25), this enhancement cancels the suppression factor Y^Yy 
that comes from the term Tb/Mb, and we have 0(1) values of ei, 62. This is the heart of resonant 
leptogenesis naturally realized with the minimal flavor violation. 

The Boltzmann equations that govern the yields of the singlet neutrinos Ni,N2 and the 
B — L charge are written as 



dY/ 



dz 

dz 

dYB-L 

dz 



sH{yvB-L) 



sH{yvB^L) 



sH{yvB-L) 



Yj 



iVi 



^ - 1 7D, + 



Yj 



N2 



^N2 



- 1 7D2 + 



El 1 Yb^l _ 



a=l 



2 Yl" 




(29) 
(30) 
(31) 
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' ^l'^ respectively denote the yields of the singlet neutrino A^^ and the lepton doublets LjS, 
and IKjjvb-l) denotes the Hubble constant at the temprature T = i/vb-l- iDa is the space- 
time density of the decay events of Na, and 7^/ the space-time density of the processes mediated 
by the U{1)b-l gauge boson: NaNa ^ //, where / indicates a standard model fermion. We 
neglect other terms that give negligibly small contributions to the yields, in accordance with 
ref. [7]. The general form of the Boltzmann equation and expressions of 7's are given in ref. 

At small z, U{1)b~l gauge interaction keeps the singlet neutrinos in thermal equilibrium. 
Below the sphaleron decoupling temperature, the sphaleron process ends and the B — L charge 
density is no longer converted to the baryon number density. Noting these, we numerically 
solved the Boltzmann equations from z = 0.001 to z = yvB-h/Tsph with the initial conditions: 

y^,(2 = 0.001) = F^^(2; = 0.001) , (32) 
Yb^l{z = 0.001) = , (33) 

where we simply assumed that there is initially no B — L charge density. The resultant B — L 
charge density is converted to the baryon number density through the sphaleron process [T6] . 
The yield of the baryon number, n^/s, is related to Yb-l by the following formula: 

— = '^yB-liz = yVB-L/Tsph) ■ (34) 



5 Numerical Analysis on Resonant Leptogenesis 

Solving the RG equations in Section 3 and the Boltzmann equations in Section 4, we calcu- 
late the baryon number density for various values of the parameters, and compare it with the 
observed value eq. (1). The parameters are 

The U{1)b~l gauge boson mass (at the scale Mq): Mz' = 2\^gB~LVB~L 
The U{1)b-l gauge coupling constant (at the scale Mq): Qb-l 
Scale above which the flavor symmetry of two singlet neutrinos appears: A 
Approximate value of the mass of the singlet neutrinos (at the scale Mq): uvb-l 
The parameter of the Yukawa kernel matrix O: a (g C) 
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We set the U{1)b~l breaking scale Mq at yvB~L, and the sphaleron decoupling temprature 
Tgph at 150 GeV [17]. We adopt the following values: 

A = 10^ GeV , 

%^ = 0.005 . 

while substituting various values into yvB-i and Mz'- We assume the normal hierarchy of the 
light neutrino masses. With these parameters, we make contour plots of the baryon number 
yield on the plane of real and complex parts of a that parametrizes the neutrino Dirac Yukawa 
matrix Yij. The baryon number density depends on a only through the form of Y^Yy^ (17). 
Since Y^Yy is invariant under the transformation: i?e[a] — i?e[a] + vr, we restrict the range 
of a to < Re[a\ < n. 

We also make contour plots of the baryon number yield for A = 10^^ GeV, for = 0.02 
and for the inverted hierarchy of the light neutrino masses, to study the impacts of these pa- 
rameters on the result. 

In Figure 2, we show how the CP asymmetry parameters, ei, €2, evaluated at the scale Mq 
depend on A. Here we set g%_i^/An = 0.005, yvs-L = 2 TeV, Mz' = 5 TeV and a = 2 + i. 
We see that above the scale of 10^ GeV, the CP asymmetry parameters of order 0.1 are au- 

0.14 

0.12 

0.10 
^ 0.08 
^ 0.06 

0.04 

0.02 

0.00 

4 6 8 10 12 14 

Logio(A/GeV) 

Figure 2: CP asymmetry parameters, ei (solid) and €2 (dashed), at the scale Mq as functions 
of the singlet neutrino symmetry scale, A. 

tomatically realized by the RG evolutions. Also, they only weakly depend on the scale of A, 
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which means that the baryon number yield is not affected by the scale above which the ffavor 
symmetry of singlet neutrinos appears. 



On the left and right of Figure 3, 4, 5 and 6, we show contour plots of the baryon number 
yield, ub/s, on the plane of (-Re[a], /m[a]) for (M^/, yvB-i) = (3, 1), (3, 2), (4, 1.5), (4, 3), (5, 2), 
(5, 3) TeV. Thick solid (red) contours denote the regions where the resultant baryon number 
yield fits with the observed value eq. (1). Cases with {Mz',yvB~L) = (3,1), (4,1.5), (5,2) 
TeV are particularly interesting because the singlet neutrinos can be produced by the decay of 
U{1)b~l gauge boson if it is produced as a resonant state at colliders. 
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Figure 3: Contour plots of the baryon number yield, tib/s, as a function of the complex 
parameter a that parametrizes the neutrino Dirac Yukawa matrix Y,y. We take Qb-l/'^'^ = 
0.005 and A = 10^ GeV, and assume the normal hierarchy of the light neutrino masses. Left: 
{Mz',yvB^L) = (3,1) TeV, Right: {Mz' ,yvB~L) = (3,2) TeV. The numbers on the contours 
denote the corresponding values of logiolns/s). The thick red line denote the region where the 
baryon number yield fits with the observed value eq. (1). 



Figure 6 is the same contour plot as the left of Figure 5, but we here take A = 10^^ GeV. 
Figure 7 is the same contour plot as the left of Figure 5, but we here take g\_j^/4:'K = 0.02. 
Figure 8 is the same contour plot as the left of Figure 5, but we here assume the inverted 
hierarchy of the light neutrino masses. 

We notice from the left panels of Figure 4 and 5 that a sufficiently large amount of the 
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Figure 4: Same as Figure 3, but we here have 
{Mz',yvB-L) = (4,3) TeV. 
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Left: {Mz',yvB-L) = (4,1.5) TeV, Right: 






1 1 ' 

^ 1 ; ■ 


'! / \ \ 

'II V ; 

I ; \ 1 

II 1 I 
1 1 1 
1 1 1 

1 \ ' 1 

I \ 





Figure 5: Same as Figure 3, but we here have Left: {Mz',yvB-L) 
iMz',yvB-L) = (5,3) TeV. 



(5,2) TeV, Right: 
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Figure 8: Same as the left of Figure 5, but we here assume the inverted hierarchy of the hght 
neutrino masses. 

baryon number yield can be produced in benchmark points with (M^/, yvs-i) = (4, 1.5) TeV 
and (5,2) TeV, and g\_j^/4:Ti = 0.005, which are of interest for collider studies. Moreover 
the observed baryon number yield, eq. (1), is obtained for moderate values of (i?e[Q;], /m[a;]), 
which implies that no fine-tuning is required for the neutrino Dirac Yukawa coupling in order 
to reproduce the observed value. 

From Figures 3, 4 and 5, we see that the heavier the singlet neutrinos are, the more baryon 
number yield we have. This is simply due to the fact that leptogenesis ends at the sphaleron 
decoupling temperature, which is assumed to be Tgph — 150 GeV. If the singlet neutrino mass is 
smaller, less singlet neutrinos decouple from the thermal bath and do out-of-equilibrium decay 
before the radiation of the Universe cools down to the sphaleron decoupling temperature. We 
also notice that we have more baryon number yield when Mz' < 2yvB~L than otherwise. 
This is because in the latter case, an on-shell U{1)b-l gauge boson mediates the processes: 
NaNa Z' ^ ff {Na- siuglct ueutriuo, /: standard model fermion), and singlet neutrinos are 
washed out more efficiently. 

From Figure 5 and 6, we confirm that the baryon number yield hardly depends on the 
scale A, above which the fiavor symmetry of the singlet neutrinos appears. This is in accord 
with Figure 2, which shows that the CP asymmetry parameters vary only by factor 2 with the 
change of A. 

Comparing Figure 7 with Figure 5, we see that the baryon number yield decreases for larger 
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Qb-Li which reflects the fact that the washing-out of singlet neutrinos is more efficient when 
the U{1)b-l gauge couphng is stronger. 

Finally from Figure 5 and 8, we see that the baryon number yield drastically decreases with 
the inverted hierarchy of light neutrino masses. This is because the neutrino Dirac Yukawa 
coupling squared, eq. (18), cannot have large CP- violating phase when two of the mass eigen- 
values for the light neutrinos are nearly degenerate, as in the inverted hierarchy case. 



6 Conclusions and Discussions 

We have studied the resonant leptogenesis scenario in the model where SM is extended 
with U{1)b~l gauge symmetry which is broken by the Higgs mechanism at the TeV scale and 
with TeV scale singlet neutrinos. We further assume a flavor symmetry on two generations 
of the singlet neutrinos which is explicitly broken by the neutrino Dirac Yukawa coupling at 
a scale above TeV. Renormalization group evolutions generate a splitting of the singlet neu- 
trino Majorana masses that is proportional to the neutrino Dirac Yukawa coupling. Hence 
the CP asymmetry parameters for the singlet neutrino decay process are automatically en- 
hanced to order 0.1 by the resonance, even though the neutrino Dirac Yukawa coupling is much 
smaller. With the above setup, we have calculated the baryon number yield by solving the 
Boltzmann equations, taking into account the effects of the U{1)b~l gauge boson that can 
wash out the singlet neutrinos. We parametrize the neutrino Dirac Yukawa coupling by the 
light neutrino mass eigenvalues obtained from the neutrino oscillation data and an arbitrary 
2x2 complex-valued orthogonal matrix. We have derived the baryon number yield for various 
values of the singlet neutrino mass and U{1)b~l gauge boson mass, and compared it with the 
observed value, eq. (1). We have shown that a sufficiently large amount of baryon number 
yield is generated with several TeV masses for the singlet neutrinos and U{1)b~l gauge boson. 
The most notable result is that the right amount of baryon number yield can be obtained for 
(M^/, yvB-L) = (4, 1.5) TeV and (5, 2) TeV with ^|_i/47r = 0.005, in spite of the fact that an 
on-shell U {1)b~l gauge boson washes out singlet neutrinos in the course of thermal leptogenesis. 

Let us discuss the implications of the above result on collider physics. At the LHC with 
y/s = 14 TeV, the cross section for a U{1)b-l gauge boson that eventually decays into 
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is given by 



a{pp ^ Z' ^ = 0.200 fb for Mz' = 4 TeV and ^|_i/47r = 0.005, (35) 

a{pp ^ Z' ^ = 0.0271 fb for Mz' = 5 TeV and qI^l/A-k = 0.005, (36) 

where we have integrated over the muon invariant mass from Mz' — 1 TeV to M^/ + 1 TeV. SM 
background is neghgibly small in these regions. We therefore expect to discover the U{1)b-l 
gauge boson with the integrated luminosity below 100 fb~^ in these cases. The production cross 
section for the two generations of singlet neutrinos is given by (z = 1, 2) 

a{pp ^ Z' ^ N,Ni) = 0.138 fb for {Mz', VVb-l) = (4, 1.5) TeV and ^|_i/47r = 0.005, 



a{pp ^ Z' ^ NiNi) = 0.0181 fb for (M^/, vvb-l) = (5, 2) TeV and gl^L/^n = 0.005. 



A singlet neutrino decays into a light neutrino and the Higgs boson or the longitudinal com- 
ponent of Z boson, or into a charged lepton and the longitudinal component of W boson. By 
combining several decay modes, it is possible to identify singlet neutrinos [TS] with the inte- 
grated luminosity of several hundred fb~^ even in the case of (M^/, i/vb-l) = (5,2) TeV. 
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Appendix 

The number density of a particle ip (a right-handed neutrino in our case) with mass 
in the early universe is evaluated by solving the Boltzmann equation of the form [13], 



(37) 



(38) 



dz 



sH{m^) 



z 




■1 



(?/; + a + ^ « + J + ■ ■ ■ ) 



YY 



7 



(z+j + ^V^ + a + -- -) 



(39) 



yeqyeq 
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where = n^/ sis the ratio of and the entropy density s, 2; = and H{m^) is the Hubble 
parameter at a temperature T = m^. The right hand side of eq. fl39|) describes the interactions 
that change number of ip, and 7^"^ is the space-time density of scatterings in thermal equilibrium. 
For a dilute gas we take into account decays, two-particle scatterings and the corresponding 
back reactions. One finds, for a decay the particle ip, 



(40) 



where Ki and K2 are the modified Bessel functions, and Trs is the decay width. For two body 
scattering one has 



'y{'ip + a^i+j-\ ) 



(41) 



where s is the squared center-of-mass energy and the reduced cross section a for the process 
■ip + a-v^i + j + -- - is related to the usual total cross section a{s) by 



(42) 



The total reduced cross section for the process f + f ^ N + N mediated by the Z' boson 
/ denotes the SM fermions) is given by 



^ . , 1047r 2 



X 



{x - yf + yc 



(x-4) 



(43) 



where y = and c = [Tz'/Mi ) with the decay width of Z' boson 



6 



3 

3 ( 1 _ 1 y _ 4) + 13 



(44) 
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